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The nebular spectrum of the type la supernova 2003hv: 
evidence for a non-standard event 
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ABSTRACT 

The optical and near-infrared late-time spectrum of the under-luminous Type la 
supernova 2003hv is analysed with a code that computes nebular emission from a 
supernova nebula. Synthetic spectra based on the classical explosion model W7 are 
unable to reproduce the large Fe Ill/Fe II ratio and the low infrared flux at ~ 1 year 
after explosion, although the optical spectrum of SN 2003hv is reproduced reasonably 
well for a supernova of luminosity intermediate between normal and subluminous 
(SN 1991bg-like) ones. A possible solution is that the inner layers of the supernova 
ejecta (v < 8000 km s _1 ) contain less mass than predicted by classical explosion models 
like W7. If this inner region contains ~ 0.5M Q of material, as opposed to ~ O.9M in 
Chandrasekhar-mass models developed within the Single Degenerate scenario, the low 
density inhibits recombination, favouring the large Fe in/Fe n ratio observed in the 
optical, and decreases the flux in the Fe II lines which dominate the IR spectrum. The 
most likely scenario may be an explosion of a sub-Chandrasekhar mass white dwarf. 
Alternatively, the violent / dynamical merger of two white dwarfs with combined mass 
exceeding the Chandrasekhar limit also shows a reduced inner density. 
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1 INTRODUCTION 

Type la Supernovae (SNela) are the most homogeneous 
type of stellar explosions, and one of the most luminous. 
They are used to constrain the energy content of the Uni- 
verse and provided the first direct evidence of acc elerated 
expansion (|Riess et al.lll998l : |Perlmutter et alJll999h . SNela 
are not standard candles per se, but the y can be stan- 
dardised (with a dispersio n of ~ 0.15 mag: I Germany et al.l 
12004 iFolatelli et aiTl2010h . thanks to a rather tight rela- 
tion between their luminosity and t he shape o f thei r light 
curve. This was first suggested by IPhiHipsl i| 1993ft , who 
parametrised the light curves of SNela based on a single 
quantity, Ami5(B), the decline of the B band magnitude 
from maximum to 15 days later. 

SNela come with a rather large range of lumi- 
nosities. Their absolute magnitudes range from ~ —17 
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for the dimmest events like SN 1991bg (|Filippenko et all 
1 19921 : iLeibundgut et al J 1 1993ft to ~ —20 for the superlu- 
minous, "sup er- Ch andrasekhar" SNe llHowell et all 120061: 
Hicken et aiTl2007l ; IScalzo et ail 1201(1 iTaubenberger et all 



20111 ). but the bulk of spectroscopically normal SNe has 
a smaller spread, between ~ —18.5 and ~ —19.5. These 
are the majority of SNe la. At the luminous end of this 
range confusion with the spectroscopically peculiar 1991T- 
like events may occur, since these show similar decline rates 
and lumino sities as some spectros copically normal SNe (e.g. 
SN1999aa, iGaravini et al.1 |2004) . At the dim end, how- 
ever, peculiar events similar to SN 1991bg, with Ami5(B) ~ 
2 mag appear to be separated from even the dimmest of 
the normal SNe, such as SNe 1992A or 2004eo, which have 
Ami5(B) ~ 1.4 mag. Only few SNe are known that hav e 
intermediate decline rates (e.g. ITaubenberger et al.| [2008). 
A classi cally known examp le is SN 1986G (Amis(B) ~ 
1.8 mag, IPhillips et al.l [1987). which however shows many 
properties of the SN 1991bg group, in particular the pres- 
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ence of a strong Ti n absorption trough between 4000 and 
4500 A at maxim um light and the ev olution of the Si II ab- 
sorption velocity (|Benetti et al.ll2005l ). 

Since 1991bg-like SNel a are almost exclu sively ob- 
served in early-type galaxies (|Hamuv et all boOO) the ques- 
tion arises whether they may have entirely different progen- 
itors. For the bulk of SNela a carbon-oxygen (CO) white 
dwarf accreting hydrogen in a binary system has tradition- 
ally been the scenario of choice, because it is thought to be 
able to give rise to homogeneous explosions. The observed 
range of SN luminosities can be attributed to the synthe- 
sis of different amounts of 56 Ni in the explosion (Arnettl 
19821; iBranchll 19921; iRiess et al.|[l996l;ICappellaro et al.lll997F 



Contardo et aLlfeood : iMazzali et al.ll2007al : IStritzinger et al] 
2006). The most promising way to tune this parameter is 



given b y the explosion scenario known as Delayed Deto- 
nation (|Khokhlovlll99ir ). although the details of the tran- 
sition from subsonic to supersonic burning are not clear. 
A recently revived alternative scenario is that some frac- 
tion of the SNela come from the merging of two white 
dwarf s whose combined ma s s exceeds the Chandrasekhar 
mass (|lben fc Tutukovlll984l ; I Webbinkl [l98i ; IPakmor et ail 
2010). Another scenario is the explosion of a white dwarf 
with a mass below the Chandrasekhar limit. This is thought 
to occur if the white dwarf accretes helium from a com- 
panion, the helium shell which builds up is very unstable 
and can detonate, triggering the explosion of the entire star 
(Wooslev fc Weaverlll994l ; iLivne fc Arnettlll995l ). This sce- 
nario was not favoured be cause it was report ed to lead 
to inconsi s tent S N spectra l|Nugent et al.lfl997l ). However, 
ISim et all (|2010T ) present sub-Chandrasekhar SNIa models 
where the overlying He layer is omitted and suggest that 
these models reproduce the observed range of light curve 
properties of SNela (from intermediate to dim), and that 
spectra obtained from their models also resemble observed 
spectra. 

Still, most SNela, excluding perhaps the least lumi- 
nous, 1991bg-lik e ones, seem to be con sistent with the Chan- 
drasekhar mass (Mazzali et al.ll2007al 1. 

A powerful method of investigating the properties of SN 
ejecta is late- time spectrosc o py, w hich probes the inner lay- 
ers of the SN. lMazzali et all (20101 ) showed that the simulta- 
neous availability of optical and infrared (IR) late-time spec- 
tra makes it possible to account for all the most abundant 
elements produced in a SN, and to infer their abundances. 
Not many late-time IR spectra are a vailable for SNela. One 
of the best examples is SN 2003hv. iMotohara et all (|2006l ) 
published IR spectra which iLeloudas et al.l (|2009h rescaled 
to match IR photomet ry taken at an epoch similar to that 
of a mid-IR spectrum l|Gerardv et al.ll2007l ). 

This is not the only thing that makes SN2003hv in- 
teresting. The SN ha s a decline rate Ami5(B) = 1.61 mag 
(Leloudas et all 120091 ). A SNIa with such a decline rate 
could be the missing link between spectroscopically nor- 
mal SNela and underluminous SNe of the 1991bg class. 
SN 2003h v may be the first wel l-observed SN with this de- 
cline rate. ILeloudas et alj (20091 ) presented a very good data 
set covering the post-maximum phase. 

Additionally, SN 2003hv was noted for the peculiar 
shifts of several nebular IR emission lines, some of which also 
seem to show flat tops (ICerardv et al.ll2007l ; [Motohar a et al.l 

Oil 



similar displacements are observed in the optical lines of 
[Fe n]. This was interpreted as possibly the indication of a 
high-density neutron-capture region, where stable 58 Ni and 
54 Fe rather than radioactive 56 Ni are produced. This region 
woul d peculiarly be displ aced from the centre of the explo- 
sion l|Maeda et alj|2010al). This findi ng led to further work 
in this direction ( Maeda et alJlioiObl ). which is however not 



the subject of the p re sent p aper. 

ILeloudas et al.l (20091 Fig. 7) showed a synthetic 
optical/nea r-IR spectrum based on the classical explosion 
model W7 l|Nomoto et all 1 19841 ). This is a Chandrasekhar- 
mass model, with kinetic energy -Ekin~ 1.3 x 10 51 erg, 
which was developed to ac count for the bulk properties of 
SNela (Branch et all 1 19851 ) and has been used as a stan- 
dard reference to describe this class of SNe. While the 
overall reproduction of the spectrum is good, they en- 
counter ed a few problems, w hich they mentioned and dis- 
cussed (|Leloudas et al.ll2009l ). Most evidently, the near-IR 
flux (A > 7500 A) is overestimated by the model by a factor 
of a few. Additionally, the ratio of the two strongest optical 
features, the Fe-dominated lines near 4700 and 5200 A, is not 
reproduced. While the bluer feature is dominated by [Fe ill] 
lines, the redder one is predominantly due t o [Fe ill transi- 
tions. The bluer feature is well reproduced in ILeloudas et al.l 
(2009), while the redder one is overestimated by about a fac- 
tor of 2. Since most of the near-IR lines whose flux is also 
overestimated are themselves [Fe n] transitions this seems to 
indicate that the ionization balance in the SN 2003hv neb- 
ula is differ ent from that predicte d by W7. Among other 
possibilities, ILeloudas et al.l (20091 ) suggest that the lack of 
near-IR flux might be attributed to a sudden onset of an 
infrared catastrophe in the densest region of clumpy ejecta. 
The IR catastrophe (IRC) is a sudden shift of the nebular 
emission to mid-IR lines, which occurs when the temper- 
ature of the SN nebu l a bec omes sufficiently low (Axelrodl 
Il980l) . ILeloudas et all (2009) argue that an IRC may have 
occurred at earlier times, and that this may be possible if the 
ejecta are clumpy. They also find that positrons are likely 
to be fully trapped. 

Given the special role of SN 2003hv as a possible link be- 
tween normal and subluminous SNela, and the significance 
of the late-time data available, we have tried to reproduce 
these data with our nebular code. This paper is organised 
as follows: in Section 2 we briefly describe our nebular code; 
in Section 3 we show and discuss synthetic spectra obtained 
adopting different scenarios; in Section 4 we discuss the time 
evolution of the nebular emission; in Section 5 our results 
are recapped and discussed. Finally, Section 6 concludes the 
paper. 



2 METHOD 



l2006l ; iMaeda et al.ll2010al ). ILeloudas et alj (20091 ) find that 



ILeloudas et al.l (20091 ) present a combined optical, near- 
IR and mid-IR spectrum of SN 2003hv in the nebular 
phase. The optical spectrum was obtained at the ESO- 
VLT 320 days after B-band maximum. The near-IR spec- 
trum was obtained at Subaru 394 days after B-band max- 
imum (Motohara et alj|2006l ). while the mid-IR spectrum 
was obtained from Sp itzer 358 days after B-band maximum 
(Gerardv et al.l [2007f ) . Since the optical and IR d ata were 
not obtained simultaneously, ILeloudas et alj (20091 ) rescaled 
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the optical spectrum and the near-IR spectrum to match the 
photometry of day 355 after maximum, the epoch of the mid- 
IR observations. Given the relatively short time elapsed be- 
tween the optical and near-IR observations, it can be hoped 
that the spectrum did not evolve much in shape. We use the 
data presented in Fig. 7 of iLeloudas et al l (|2009l ). 

At these late times, the spectrum of SN 2003hv is dom- 
inated by nebular emission. The gas is optically thin to op- 
tical and IR radiation, and it is heated by the deposition 
of the 7 rays and positrons from the decay of 56 Co to 56 Fe 
( 56 Ni has almost completely decayed to 56 Co). Collisional 
heating is balanced by cooling via emission, and the energy 
is released in a number of line transitions, most of which are 
forbidden. 

We use a non-local thermodynamic equili brium (NLTE) 



by 



nebul ar code which was developed 
l|200lh_ and is based on ideas 

and Ruiz-Lapuente fc Lucvl lll992 | ). Although a t hree- 
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dimensional version is available jlVIaurer et al.l l2010al ) the 
code was used in its one-dimensional version here. Our (one- 
dimensional) models consist of about 20 radial shells of vary- 
ing density and composition. The emission and deposition of 
7 rays an d positrons is computed with a Montecar lo scheme, 
following lCappellaro et all (I1997T ) and iLucvl d 2005). using ef- 
fective opacities for 7 rays and positrons. 

After energy deposition, an iteration process is initi- 
ated. Starting with an initial guess for the ionisation state 
and the temperature of the gas, non-thermal ionisation and 
excitation rates are derived using the Bethe and the optical 
approximation respectively jAxel rod 1980; iRozsnvai et aJj 
1980; Ma urer fc MazzalU l2010h . which have b een shown to 
be acc urate to at least 20% for H and He (|Maurer et al.l 
2010b). For other elements the accuracy of these approxima- 
tions is difficult to estimate, since there is a strong depen- 
dence on poorly known atomic data. Since our nebular code 
neglects radiation transport, there is no photo-ionisation. 
The collisional ionisation rates are balanced with the radia- 
tive and the di-electronic rec ombination rates taken from 
lAldrovandi fc Peouignotl (|l973T ). Charge exchange reactions 
are not considered. After obtaining the new ionisation state 
and consequently a new electron density, the level popula- 
tion of each ion is derived by solving a matrix containing 
collisional (de-)excitation and radiative de-excitation rates. 
The radiative rates are reduced according to the Sobolev 
optical depths of the corresponding lines. Knowing all the 
level populations, a radiation field can be generated and can 
be compared to the total deposited luminosity. The electron 
temperature (and consequently the electron density) is then 
adjusted in each iteration step, until the total luminosity of 
the radiation field has converged to the deposited luminos- 
ity. 

The code can treat the neutral and singly ionized states 
of H, He, C, N, O, Ne, Na, Mg, Si, S, Ar & Ca and the 
singly and doubly ionized states of Fe, Co & Ni. For light 
and intermediate-mass elements typically the lowest 10 — 20 
levels are taken into account in the excitation calculation. 
For iron-group elements roughly the 100 lowest energy levels 
of each ion are considered. Collisional data are available for 
several hundred transitions of H, He, C, N, O, Ne, Na, Mg, 
Si, S, Ar, C a, Fe, Co fc Ni . Most of the collisional data are 
taken from lAxelrodl l|l980T ) or from online databases (e.g. 



TIPTOFfl A. PradharQ) and are updated frequently. Col- 
lisional data are available for the most important forbidden 
lines, but those for Co are unfortunately very poor. If not 
available, the c ollision strengt h of forbidden lines is crudely 
set to O.Olgrigfu (|Axelrodl ll980). where g\ and g u are the sta- 
tistical weights of the lower and the upper atomic level re- 
spectively. Unknow n collisional data for al lowed transitions 
are obtained in the Ivan Regemorter] l|l962T l approximation. 

Although missing collisional data pose a problem for all 
kinds of nebular calculations, the largest uncertainties in our 
calculations are caused by neglecting photo-ionization. It is 
clear that the neutral species of the iron-group elements are 
photo-ionized almost completely in SNe la at the epochs of 
interest (200 - 400 days after the explosion). To simulate this 
effect, Fe I, Co I and Ni I are not included in our calculation. 
The ratio of Fe ill to Fe 11 is determined by the non-thermal 
electron ionisation rates, which is a good approximation. For 
light and intermediate-mass elements the problem is more 
serious and it is not clear if we can produce the ionisation 
state of those elements correctly. However, since our syn- 
thetic spectra are dominated by Fe 11 and Fe ill emission 
lines, we expect them to be reliable within the uncertainties 
of the atomic data. 

The code has been used to investigate the morphol- 
ogy and composition of t he spectra of various SNe (e.g. 
iMazzali et all2005l . l2007bh . A 1-zone version makes it possi- 
ble to derive very quickly the basic properties of the emitting 
nebula (mass, energy, composition) using the maximum ob- 
served emission line velocity as an outer boundary. In a more 
sophisticated stratified version the density varies as a func- 
tion of velocity (which is equivalent to radius in the homol- 
ogously expanding SN nebula) and so do the abundances. 
These can be modified in order to reproduce the observa- 
tions, yielding a rather accurate description of the structure 
of the SN ejecta. The densit y structure can be t aken from a 
model calculation (e.g. W7, Tanaka et al.ll201ll ). but it can 
also be derived from the emitted luminosity and the line 
shapes. 

For all calculation s we adopt for SN 200 3hv a distance 
modulus of 31.58 mag (jLelou das et al. 2009) and a redden- 
ing E(B -V) = 0.016 mag dSchlegel et al.lll998l). Reddeni ng 
within the host galaxy is negligible I Leloudas et al.|[2009T l. 



3 MODELS 

The late-time optical spectra of SNe la are dominated by 
strong emission in forbidden lines at ~ 4700 and 5200 A. 
The features are broadened because they are blends of many 
lines, so that their width is larger than the real velocity of the 
emitting gas. The bluer feature is dominated by [Fe ill] lines 
(AA4607, 4658, 4667, 4701, 4734, 4755, 4881), with some 
[Fc 11] lines (AA4814, 4874, 4890) also contributing, while the 
redder one has comparable contributions of [Fe 11] (AA5111, 
5159, 5222, 5262) and [Fe ill] lines (AA5270, 5412). A cor- 
relation exists between the width of the features and the 
luminosity of the SN as described by the light-curve decline 
paramater Am\${B) l|Mazzali et al.lll998h . This is because a 
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Figure 1. Late-time optical spectra of 4 rapidly declining SNe la: 
SN 2003hv (Ami5(B)= 1.61 mag, spectroscopically normal) - 
black; SN2004eo (Ami^(B) = 1.47 mag, spectroscopically nor- 
mal) - green; SN 1989B (Ami5(B) = 1.28 mag, spectroscopically 
normal) - blue; SN 1986G (Ami5(_B)= 1.73 mag, spectroscopi- 
cally peculiar) - red. The spectra have been normalised at the 
peak of the 4700 A emission. 



more luminous SN burns more of the progenitor white dwarf 
into 56 Ni, which is observed at late times as Fe emission. 
The availability of the near-IR spectrum makes it possible 
to study all elements that are supposed to be abundantly 
synthesised in a thermonuclear explosion. 

To begin with, we compare the optical spectrum 
of SN 2003hv to those of other SNe la with similar de- 
cline rates. Figure 1 shows the blue part of the nebu- 
lar spectra of 4 SNe la: SN 2003h v (Amis(S)= 1.61 mag), 
SN1986G (IPhillips et all Il987l, Am 15 (B) = 1.73 mag), 
SN2004eo ilPastorello et al l 2007] . Am 15 (B)= 1.47 mag), 
and SN 1989B dwells et al.1 1 19941 . Ami 5 (B)= 1.28 mag). 



With the exception of SN 2003hv all spectra are from the 
Asiago database. They have all been corrected for redden- 
ing, and are plotted after normalizing at the peak of the 
4700 A line. SN 1986G has the narrowest emission lines, as 
may be expected given that this SN has the most rapidly de- 
clining light curve. The other 3 SNe have all rather similar 
line width. SN 2003hv stands out for the much higher ratio 
of the 4700 A emission feature with respect to the 5200 A 
one. 

We begin with modelling the optical spectrum alone. 
For this purpose we use the optical spectrum calibrated to 
the original photometry of the epoch when it was actually 
observed, 320 days after maximum, thus avoiding any un- 
foreseen evolution between then and the epoch of the near- 
or mid-IR spectrum. We discuss here various models. 



3.1 One-zone models 

We use the one-zone version of the code in order to de- 
termine the velocity width of the emitting nebula, assum- 
ing that it is spherically symmetric. The emission lines of 
SN 2003hv can be reproduced for a width of the emitting 
region of 8000 km s _1 . We also measure the FWHM of the 
4700A line, finding a value of 13500 ± 450 km s" 1 . Both 
of the se values are consistent with Fig. 2 of iMazzali et al.l 
(|l998h . indicating that SN2003hv is located at the edge 
of the distribution of normal SNe la, close to SNe 1986G 
and 1993L. The value of Amis(_B) that would be obtain ed 
from the fit shown in Figure 2 of IMazzali et al.l l| 19981 ) is 
~ 1.4 ma g. The photometric value (Atois(-B) = 1.61 ± 
0.02 mag; ILeloudas et alj 120091 ) is different but not incon- 
sistent, given the la rge dispersion of the relation shown in 
IMazzali et ail (|l998h . 

However, SN 2003hv stands out because it shows a 
much larger ratio of the bluer, Fe ill-dominated line near 
4700 A over the redder, Fe n-dominated one near 5200 A ( sec 
Fig.l). Reddening towards SN2003hv is small. SNe 1989B 
and 1986G a re both highly reddened (E(B - V) = 0.37 
and 0.60 mag; IWells et aLllilffll ; IPhillips et al.lll999l . respec- 
tively), but even when this is corrected for a clear differ- 
ence remains. The same is true if we compare SN 2003hv 
with SN 2 004eo, which has a low er reddening (E(B — V) = 
0.11 mag. iPastorello et al.ll2007n and similar line width. 

A reasonable 1-zone model for the optical spectrum of 
SN2003hv is obtained for a 5e Ni mass of 0.42M Q and is 
shown in Figure 2. This is actually rather a large value 
for the S N decline rate : for SN2004eo, with Ami 5 (B) = 
1.47 mag. [Mazzali et al.1 (|2008D found M( 56 Ni) = O.35M . 
However, this value is consistent with the estimate of 
ILeloudas et al.l (I2009T ) based on the luminosity at peak. The 
main reason for the large 56 Ni mass is indeed the need to 
fit the strong Fe in emission. The large Fe Ill/Fe II ratio 
(~ 1 : 1) which characterizes the model is obtained assum- 
ing that no stable 54 Fe is present in the centre. This isotope 
acts only as a coolant and its effect is to reduce the ioniza- 
tion of Fe. Not introducing stable Fe keeps the ionization 
ratio high. Co and Ni have a similar ionization state as Fe. 
Intermediate-mass elements (IME) are more highly ionized. 
The mass used in the model within 8000 km s _1 is only 
0.46Mq, which is less than in 1-zone models for the other 
SNe (e.g., SN 2004eo requires ~ 0.55M Q within 7300 km s _1 ) 
and significantly les s than the prediction of W7 (0.63Mq, 
llwamoto et al.Hl999T ). So, in the 1-zone model, SN2003hv 
has a small mass at low velocities, and most of it seems to 
have been burned to 56 Ni. 

Although the 1-zone model does a reasonable job of 
reproducing the optical spectrum, it can be seen to fail al- 
ready in the reddest part: the [Fe n] emission lines redwards 
of ~ 7500 A are overestimated. This becomes even worse in 
the near-IR, as is shown in Figure 2. 

There are therefore several different lines of evidence 
indicating that SN 2003hv may not conform with the norm of 
SNe la. This can be further studied using stratified models. 



3.2 Stratified models: W7 

In the stratified version, our nebular code can be used to test 
specific explosion models, or, as we do here, to construct a 
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Figure 2. The late-time optical/near-IR spectrum of SN 2003hv compared to various synthetic spectra: a 1-zone model (red); a stratified 
model based on the W7 density distribution but with modified abundances (green); a stratified model also based on W7 but with the 
additional inclusion of clumping in density (blue). 



best-fitt ing density-abundanc e distribution by fitting obser- 
vations. [MazzafLeTlil] (|2007al ) fitted the late-time spectra of 
23 SNe la of different luminosities using the density distri- 
bution of W7. They changed the 56 Ni content by modifying 
the outer extent of this species, as guided by the width of 
the observed nebular lines of Fe. They found that not only 
the mass and distribution of 56 Ni correlates with the lumi- 
nosity: evidence that the distribution of IME such as silicon 
and sulphur is complementary to that of nuclear statistical 
equilibrium (NSE) material suggests that most of the ma- 
terial which is not burned to Fe-group elements must have 
been burned to IME. This supports the notion of a common 
mass for all SNe la, or at le ast for most of them, excluding 
perhaps 1991bg-like events (|Hachinger et al]|2009h . 

Here we test whether the W7 density distribution can be 



comp atible with SN 2003hv. The W7 model (|Nomoto et all 
1 1984 ) was created to fit the basic properties of average 
SNe la. The underlying assumption is that SNe la result from 
the deflagration of Chandrasekhar-mass C-O white dwarfs. 
We now know that this is at best a second ary channel to 
the production of SNe la (|Sahu et al.| [2008). However, the 
density- velocity distribution of W7 is representative of typ- 
ical SN la ex plosions: delayed deto nation models are indeed 
very similar (|lwamoto et al.lll999l ). 

Using the W 7 density struc t ure, an d modifying only the 
abundances as in lMazzali et al.l ((200 7a). it is possible to get 
a reasonable fit of the optical spectrum of SN 2003hv (Figure 
2). This model has M( 56 Ni) = 0.42M Q , similar to the 1-zonc 
model. 56 Ni is located between 2000 and 10000 km s -1 , with 
abundances peaking at ~ 70 % at 5000-7000 km s" 1 (Figure 
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Figure 3. Abundances for the model based on the W7 density 
distribution. 



3). Stable 54 Fe is located only in the innermost region, be- 
low 4000 km s" 1 : the tota l mass of 54 Fe is O.O9M , as in 
W7 l|lwamoto et al.l Il999l ) . IME must already be present 
as deep as ~ 4000 km s _1 and become dominant above 
~ 8000 km s _1 , otherwise the mass of 56 Ni would be too 
large and the Fe lines too broad. Regions above 10000 km s _1 
are not sensitive to the nebular treatment as their density 
is too low. 

The synthetic spectr um basically reproduces the results 
of iLeloudas et all 12009), showing that the two codes give 
consistent results. While the blue part of the spectrum re- 
produces SN 2003hv, the emission lines in the red part of the 
optical spectrum as well as in the near-IR are obviously too 
strong (Figure 2). Most of the lines in those regions are due 
to [Fe n]. The strongest features in the bluer part of the op- 
tical spectrum, on the other hand, are either dominated by 
[Fe in] (the line near 4700 A) or have a non- negligible [Fe in] 
contribution (the line near 5200 A) . This suggests that the 
ionization degree of Fe is too low in the model. This can- 
not be improved as long as the high central density that 
characterizes W7 is kept. Another reason for the low ioniza- 
tion is that in the inner part of the ejecta stable Fe-group 
isotopes dominate. The role of these isotopes is only that 
of a coolant. Replacing stable Fe-group isotopes with 56 Ni 
would increase the ionization because of the added heat- 
ing, but it would also lead to an increase of the luminosity. 
Replacing them with IME would further increase the [Fe n]- 
[Si i] line near 1.6/j.m, which is already too strong (Figure 
2), and replacing them with ox ygen would cause a n emission 
line due to [O i] 6300, 6363 A jKozma et alj|2005t ). which is 
not observed. In any case, none of these changes would be 
justified by explosion physics, since some neutron-rich ma- 
terial should be synthesised at the high central density of a 
Chandrasekhar-mass white dwarf. 



3.3 Stratified models: clumped W7 

ILeloudas et~aH (|2009T l suggest that if the inner ejecta were 
clumped the lack of near-IR radiation may be explained be- 
cause cooling would be less efficient than in the absence of 
clumping. Our code allows the simulation of clumping in 
density: if we define a filling factor / < 1, the gas is as- 



sumed to be concentrated in clumps which fill a fraction / 
of the volume. We tested a density clumping scenario with 
a moderate filling factor of 0.5, which does not affect the j- 
ray deposition. The result of this is that the higher density 
in the clumps actually favours recombination, and further 
decreases the Fe Ill/Fe II ratio. A model with clumping is 
shown in Figure 2. The synthetic spectrum is still charac- 
terised by a high ratio of near-IR to optical flux. Therefore, 
moderate clumping in density is unlikely to explain the ob- 
servations of SN 2003hv. 



3.4 Stratified models: a modified density 

Since there seems to be no solution for the spectrum of 
SN 2003hv using a W7 density, we next investigate whether 
different density structures might lead to better results. The 
advantage of modelling nebular-phase spectra of SNela is 
that photoionization can be, to a first approximation, ne- 
glected in the iron core (see Section 2). Combined with the 
homologous expansion of the SN ejecta, this makes nebular 
spectroscopy a very powerful tool. If the properties (den- 
sity and abundances) of a certain velocity shell are modified 
while the rest of the nebula is left unaffected it is possible 
to derive the properties of that specific velocity shell. Thus 
we can test in a rather uncomplicated way whether modi- 
fications in the density structure can lead to an ionization 
ratio which favours Fe ill more than W7 does. 

Qualitatively, in order to increase the ionization of the 
gas the recombination rate should be reduced, which could 
be achieved with lower densities. Alternatively, the ioniza- 
tion rate should be increased, which can be obtained by 
reducing the ratio of stable versus radioactive NSE ele- 
ments. This suggests that the density in the regions where 
the lines are emitte d may be smaller than in W7 (see also 
iMaeda et al.l l2010al ). and that this may best be achieved 
by reducing the content of stable NSE material, as already 
found in the one-zone model (Section 3.1). Therefore we 
arbitrarily decreased the W7 density in the inner part of 
the ejecta, where the emission lines are produced. This was 
achieved mostly by removing stable NSE material. The outer 
layers, which cannot be probed in the nebular phase because 
their density is too low, were left as in W7. At every step, 
we also modified the abundances in order to obtain as good 
a match to the data as possible, both in the optical and, in 
particular, in the near-IR. We modelled the optical/near-IR 
spectrum, now rescale d to a common epoch of day 355 after 
B maximum following ILeloudas et al.l {2009 ). The modified 
model has only ~ 0.5M© of material with v < 10000 km s _1 . 
In W7, this is w 0.9M Q . 

The resulting synthetic spectrum is shown in Figures 
4 and 5 (red line). As expected, we indeed find that when 
a reduced density is used the ionization degree is increased 
with respect to the W7-based model (Fig. 4, green line). This 
leads to a reduced flux in the red and the near-IR, which 
is the main result we tried to achieve. At the same time, 
the increased ionization causes the [Fe ill] 4700 A line to 
become somewhat stronger with respect to the 5200 A line, 
where both [Fe ill] and [Fe n] lines contribute. The ratio of 
these two Fe lines is not perfectly reproduced, as we tried to 
obtain a synthetic spectrum that reproduced the observed 
one as well as possible over a wide wavelength range. A 
major deficiency is the line near 4300 A which is much too 
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Figure 4. The late-time optical/near-IR spectrum of SN 2003hv, scaled to the photometry of day 355 after B maximum (black) compared 
to various synthetic spectra. Green: a stratified model based on the W7 density distribution but with modified abundances; red: a stratified 
model with a reduced density in the inner layers; blue: a stratified model with reduced density in the inner layers but an increased density 
in the outer layers to preserve the Chandrasekhar mass. 



weak. A synthetic [Fe n] line at 4416A is present, but it is 
not very strong. [Fe n] 6555 A is also too weak. Increasing 
the strength of those li nes would lead to the overestimate of 
other [Fe n] lines, as in lLeloudas et al. (2009). The emission 
line near 1.6 /im does not sho w a flat top: this woul d require a 
three-dimensional study (see lMaeda et al1l2010ah since this 
and other features char acterized by a flat top are peculiarly 
blueshifted in the data ijMotohara et al.ll2006l ) . In our model 
7 rays deposit in the central zone, and flat tops are not 
expected. Figure 5 shows a blow-up of the red/near-IR part 
of the spectrum. 

The various density profiles tested in the paper are 
shown in Figure 6. The 56 Ni mass included in the reduced- 
density model is only O.I8M0, and the mass of stable NSE 



elements (mostly 54 Fe) is 0.05Mq, in both cases much less 
than in W7 (which has a 56 Ni mass of ~ 0.63M p) and a stable 
NSE mass of ~ 0.2MM. Ilwamoto et al.lll999h . Some reduc- 
tion of the 56 Ni mass is always seen when going from 1-zone 
models to stratified ones, and it results from the increased 
energy dep osition in the densest inner layers of a stratified 
model (e.g. iMazzali et al]|2007bh . In the case of SN 2003hv, 
however, this reduction is dramatic. It is explained by the 
significantly reduced near-IR luminosity in the new model. 
There is now a large discrepancy with the estimate from 
the maximum of the light curve, a result also obtained by 
iLeloudas et al.l (120091 ) when comparing the peak and the tail 
of the light curve. Figure 7 shows the abundance distribu- 
tion for this model, and Figure 8 the ionization of some key 
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Figure 5. Same as above, but focussing on the low-level red and 
near-IR, part of the spectrum. The two models shown are the 
stratified model with a reduced density in the inner layers (red); 
the stratified model with reduced density in the inner layers but 
an increased density in the outer layers to preserve the Chan- 
drasekhar mass (blue). 

elements, compared to the model based on the W7 density 
distribution. 

Finally, we computed a spectrum for a model where the 
mass which is missing at low velocities, compared to the 
W7 model, is located at high velocities. The mass which 
has been removed from the densest inner layers is added at 
higher velocity in this model (Figure 6) so that the width of 
the emission lines is not increased. This model produces a 
similar nebular spectrum as the one with reduced mass (Fig- 
ures 4 and 5). Nebular models cannot explore the properties 
of the outer layers. This requires m odelling the early-time 
spectra (e.g. lHachinger et alj|2009t ). Only then may a com- 
plete model for the SN be defined. 



4 EVOLUTION OF THE OPTIC AL/NEAR-IR 
SPECTRUM 

We have shown that the optical/near-IR spectrum of the 
underluminous but spectroscopically normal SNIa 2003hv 
at ~ 1 year after explosion can be better reproduced if the 
emitting mass, which is located mostly at v < 10000 km s , 
is significantly reduced with respect to what is predicted by 
a standard explosion model such as W7. Delayed detonation 
models have a similar density distribution as W 7 in the inner 
layers l|Hoflich e t al. 19951; Hwamoto et al. HlQQSl . Fig. 3). The 
reduced mass in our model is required by the relatively high 
ionization of the gas, as signalled by the unusual strength of 
the [Fe ill] optical emission lines relative to the [Fe n] lines 
as well as by the unexpectedly low near-IR flux, which again 
consists mostly of [Fe n] lines. 



T 



W7 model 




velocity [km s^'] 

Figure 6. The density structures used for the calculations: the 
W7 density distribution (green); a modified model with a reduced 
density in the inner layers (red); a model with reduced density 
in the inner layers but an increased density in the outer layers to 
preserve the Chandrasekhar mass (blue). 
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Figure 7. Abundances for the model with the modified density 
distribution. We show the part of the ejecta which can be probed 
by nebular spectroscopy. 

iLeloudas et al . (2009, Fig. 8) have shown that the con- 
tribution of the near-IR flux to the emitted spectrum in- 
creases from 5% of the UVOIR flux (|Leloudas et al.ll2009l . 
the mid-IR flux is not included in this estimate) at ~ 100- 
200 days to ~ 10 % at about 1 year, and reaches a maxi- 
mum of ~ 40 % at ~ 600 days. We have tested whether our 
model can reproduce this trend. We have computed spec- 
tra at different epochs and integrated the flux to determine 
the relative near-IR contribution as a function of time. The 
results are show in Figure 9: although our model with a mod- 
ified inner density does not reproduce the observed values 
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Figure 8. The ionization structures for the W7 density distribu- 
tion (dotted lines) and the modified model with a reduced density 
in the inner layers (fully drawn lines). In the various panels, blue 
lines represent the more highly ionised of the two species treated. 
Truncation of the lines reflects the lack of the corresponding ele- 
ment in certain velocity shells. 



exactly, it is much closer to the observed trend than the W7- 
based model. This adds support to the idea that SN 2003hv 
had a mass at low velocities significa ntly smaller than wh at 
Chandrasekhar-mass models predict (jMaeda et a l. 2010a). 

We have not extended the calculations to beyond 600 
days, because at those epochs emission is expected to shift to 
the mid-IR, and many of the relevant lines are not included 
in our code. 



5 DISCUSSION: MASS AND RISETIME 

Our best-fit model has a mass of M c j — O.5M0 inside 
10000km s _1 . The 56 Ni mass is w O.18M . By contrast, 
a model like W7 has a mass of 0.9Mq inside the same veloc- 
ity. Models based on W7, with varying abundance distribu- 
tions (in particular different 56 Ni masses) have successfully 
been used with our code to reproduce t he nebular spectra 
of SNela of different luminosities (e.g. IStehle et a.1.1 120051 : 
iMazzali et ai1l2008l : iTanaka et al.ll201lh . The difference we 
see in this case suggests that SN 2003hv may have been the 
result of a different type of explosion, with a mass that is 
distributed differently in velocity and may be different from 
the Chandrasekhar mass. 



5.1 Mass estimates and uncertainties 

Here we address some uncertainties inherent to our mass es- 
timates. The nebular model presented above has a smaller 
mass than W7 in the inner part [v < 10000 km s _1 ). It pre- 
serves the density- velocity structure of W7 at higher veloci- 
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Figure 9. The near-IR contribution to the UVOIR flux as a 
function of time: th e points are the measured fractions (from 
iLeloudas et al.l l2009). while the orange line is the result of our 
model with reduces mass. We also show the expected behaviour 
of the modified W7 model (with a reduced 56 Ni mass) which we 
use to model the spectrum of SN 2003hv. 



ties, but those layers are not probed by nebular spectroscopy. 
Our model probably somewhat underestimates both the 
56 Ni mass and the mass at low velocities, for two reasons. 
First, the low-level flux between the strong optical emission 
lines is not well reproduced: our synthetic flux in the optical 
region is lower than the observed one by ~ 20%. Secondly, 
the flux emerging in the mid-IR contributes ss (34 ± 17) % 
to the total SN lum inosity 358 days after the explosion 
l|Leloudas et al.ll2009( ). while this fraction is only ~ 3% in 
our synthetic spectrum. 

The missing continuum in our synthetic model is the 
consequence of the absence of UV radiation transport in our 
nebular code. It is known that UV ex citation and fluores - 
cence are reponsible for such continua dLi fc McCravlir996l ). 
The missing mid-IR flux is probably caused by atomic tran- 
sitions at energies of about 0.1 eV and less, which are re- 
sponsible for mid-IR emission but are partially missing in 
the nebular code. Since an increase in flux can be obtained 
by increasing either the mass of 56 Ni or the total emitting 
mass (in both cases the energy deposited in the nebula in- 
creases), both the 56 Ni mass and the emitting mass (within 
10000 km s _1 ) could be increased by a factor 1.25 ± 0.10. 
Since our model has an emitting mass of O.5M0 within 
10000 km s , when this correction is applied we get an emit- 
ting mass of O.63±O.O5M0 . As for the 56 Ni mass, in our mod- 
els we obtained 0.18M Q . This then becomes 0.22 ± O.O2M . 
Both values are still significantly smaller than the corre- 
sponding values in W7, which has a mass of O.9M0 inside 
of 10000km s _1 and a 5B Ni mass of ~ O.63M . 

Other sources of uncertainty can affect our result. 
The distance we adopted (fi = 31.58 mag) was derived 
matching the BV RI light curve of S N2003hv to tem- 
plate light curves (|Leloudas et al.l I2009T ). This is 0.2 mag 
larger than the corrected SBF distance to the host galax; 



of SN2003hy, NG C 1201 (31.37 mag. iTonrv et all 120011 
I Jensen et al.|[2003l ). If we had used the SBF distance, both 
the 56 Ni mass and the total mass would be smaller than what 
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we have derived by ~ 10% (alternatively, either the 56 Ni 
mass or the total mass might be smaller by up to ~ 20%). 
The extinction to SN 2003hv is quite small, so that our re- 
sults are basically insensitive to its exact value. 

We have tested whether the mass which is missing in 
the inner layers may be located at high velocities by comput- 
ing a model where we increased the density of shells above 
10000 km s _1 so as to preserve the Chandrasekhar mass of 
the ejecta. The corresponding spectrum, which is shown in 
Figure 4, is not easy to distinguish from that of the model 
with a lower mass. This is because, even with an increased 
mass, the density in the outer layers is still too low to trap 7 
rays and produce line emission at the epoch of the nebular 
spectrum of SN 2003hv. 

An abundance tomo graphy experiment su ch as that 
performed for SN 2005bl (fHachinger et all 120091 ) would be 
required to determine the mass and energy in the outer lay- 
ers. We plan to perform such an analysis in the near future. 



5.2 



3 Ni mass and risetime 



iLeloudas et all (|2009h derive a 56 Ni mass o f ~ O.4M fr om 
the light curv e peak applying Arnett' s rule l|Arnettll 19821 ) as 
calibrated bv lStritzinger et al.l (|2006h using the formula 



L = 2 x 10 43 M( 56 Ni)/M 



[erg ; 



(1) 



which is calibrated for a risetime t r = 19 days. They per- 
form a simple light curve analysis and find that this value 
is however inconsistent with the late part of the light curve, 
which requires a 56 Ni mass of ~ O.2M0. They also show a 
synthetic nebular spectrum which is much brighter than the 

obs ervations. 

ILeloudas et~ai1 (2009) suggest various possibilities to 
explain this discrepancy, including an early near-IR catas- 
trophe which causes the flux to shift to the far-IR. An early 
IR catastrophe seems not to be supported by the data, be- 
cause we have seen that the gas is still quite highly ionised 
even at later times, and is charachterised by temperatures 
well in excess of 1000 K. An IRC occurring only in the in- 
nermo st layers, a possibility suggested by iLeloudas et al.l 
(2009), would probably not transfer enough energy to the 
far-IR. 

On the other hand, our model with a reduced mass at 
low velocity is capable of r eproducing the late-tim e flux for 
a 56 Ni mass similar to what ILeloudas et al. I (|2009h find from 
the late light curve. This leaves open the question of why 
such different 56 Ni masses can be derived from the early 
and the late data, wh ich usually are in good agreement 
jStritzinger et a.f]|2006h ■ Although Arnett's rule may not be 
highly accurate, a factor of 2 difference is a large discrep- 
ancy. 

We suggest that a simple solution to the problem can 
be found if the density in the inner part of the ejecta is re- 
duced. The 56 Ni mass was estimated from the peak of the 
light curve using the formula above, which assumes that 
light curve maximum occurs 19 days after the explosion 
(|Stritzinger et ai1l2006l ). However, if the light curve risetime 
of SN 2003hv was significantly smaller than 19 days, this 
discrepancy may be resolved. Using the dependence of the 
lu minosity on the 56 Ni mass and the risetime as given in Eq.5 
of IStritzinger fc Leibundgutl (|2005l ) yields the curve shown 
in Figure 10. Adopting the estimate of the 56 Ni mass from 




Figure 10. Effect on the risetime of a modified estimate of the 
56 Ni mass, including both 56 Ni and 56 Co decay. 



th e nebular model s (0.22 M(T)) instead of the value obtained 
by ILeloudas et al.l ((2009) from the peak of the light curve 
and an estimated risetime of 19 days (0.4M Q ), i.e. reducing 
the 56 Ni mass by (0.22/0.40) = 0.55, a risetime t r » 10 days 
would be inferred. This would reflect the rapid decline of the 
light curve, while the 56 Ni mass estimate from the nebular 
spectrum would hardly be affected by this choice of rise- 
time, since the decay time of 56 Co, which powers the SN 
light curve at these advanced stages, is much longer. 

A risetime of ~ 10 days is short for SNela. Estimates 
of the typical risetimes in the literature have decreased from 
« 19.5 days, a resu lt based on extrapolating the rising par t 
of the light curve dRiess et al.lll999f; lAldering et alj|2000h . 
to values bet ween ~ 19.1 da y s dConlev et al.l I2006T) and 
» 17.5 days (|Garg et al.ll2007l ; lHavden et alj|2010| ). based 
on direct observations of SNe. Both of these values are sig- 
nificantly larger than what would be required in order to 
bring the early- and late-time estimates of the 56 Ni mass 
in SN 2003hv into agreement. On the other hand, methods 
based on fitting the rising part of the light curve yield shorter 
risetimes, whic h are possibly correlated to the decline rate 
(|Strovinkll2007h . when only SNe with a sufficiently well sam- 
pled pre-m aximum photometry are an alysed. The fas test 
decliner in lHavden etal] (|201Ch and IStrovinkl (I2007T) is 
SN2004eo (Ami 5 (B)= 1.47 mag; |Pastorello et al ] |2007l ). for 
which they estimate t r — 16.6 days. However, shorter rise- 
times have been inferred or observed. SN 1994D (A mis(_B) = 
1.34 mag) is quoted to h ave t r — 15. 4 days l|Strovinkl 
l2007l ; lHavden et ail |2010T >. lHavden et all (|2010l ) observed 
a few SNe with t r ~ 13 — 15 days. The incidence of 
rapid risers is much higher for dimmer SNe (Amis^B) > 
1.5 mag), although these SNe are much more seldom seen 
than brighter o nes. Extrapolating the relation in Figure 
2 of IStrovinkl l|2007h to the decline rate of SN 2003hv, 
Ami5(B) — 1.61 mag, a risetime of ~ 13 — 15 day s is ex - 
pected, in agreement with the result of lHavden et al.l (|2010t ). 

A further suggestion of a short risetime of SN 2003hv 
comes from an attempt to compute the light curve proper- 
ties in an approximate way as done in lMazzali et al.l (|2007al ) . 
Adopting a total mass of 1.12M0, and assuming that this is 
distributed for 0.22M Q as 56 Ni, 0.07M Q as stable Fe-group 
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Figure 11. Two synthetic bolomctric lig ht curves compared t o 
the UVOIR light curve of SN 2003hv (dots, iLeloudas et al.ll2009t) . 
The fully drawn line is a model computed for the density struc- 
ture and 56 Ni mass derived from the nebular spectrum models. 
It has M= I.QMq, M( 56 Ni)= 0.18M Q , and E kin = 1.24 10 51 erg. 
The dashed line is a model where both the 56 Ni mass and the 
total mass located at velocities v < 8000 km s _1 are multi- 
plied by a factor 1.25, reflecting inaccuracies of the nebular- 
epoch modelling (see Section 5.1). This model has M= 1.12Mq, 
M( 56 Ni)= O.22M , and E Un = 1.26 10 51 erg. Both models have 
risetimes of ~ 14 days. 



material, and that most of the remaining mass located in- 
side of 11000 km s~ 1 (O.4M ) has been burned to IME, this 
yields a nuclear energy release of fa 10 51 erg. Subtracting 
from this the binding energy of a white dwarf of ~ 1.12M , 
~ 2 ■ 10 50 erg, an explosion energy -Ekin ~ 8 • 10 50 erg is 
obtained. Combined with an es timate of the opacity as dis- 
cussed in lMazzali et al.l (|2007al ). this yields a characteristic 
bolometric LC width of ~ 15 days. Looking at Fig. 2 of 
iMazzali et al.l l|2007al ). this would correspond to a peak lu- 
minosity of only ~ 3T0 42 erg s _1 , which is much less than the 
observed value ~ 10 43 erg s _1 . Therefore SN 2003hv does not 
seem to comply with the correlations established for SNe la: 
the 56 Ni mass derived assuming a relation to the peak lumi- 
nosity calibrated for a risetime of 19 days is too high for the 
decline rate of the SN. 

We finally tested the possible risetime for SN2003hv 
using the same grey Montecarlo light curve code which 
could successfully re produced the bolometric light curves 
of other SNe la (e.g . IStehle et ai]|2005l ; IMazzali et al.ll2008l ; 
iTanaka et all 1201 if ). We use the "reduced density" model 
derived from nebular spectroscopy and shown in Fig. 6. The 
model has a mass of O.5M inside 10000 km s" 1 . Outside 
of this velocity we simply assume the W7 density distribu- 
tion. The total mass is therefore 1.12M , the 56 Ni mass is 
0.22M Q , and the kinetic energy E kin = 1.26 ■ 10 51 erg. We 
also used the sam e definition of the op acity as a function of 
chemical species (jMazzali et al. I l2007al ) as in those papers. 
The resulting bolometric light curve is shown in Figure 11. 
The risetime of the synthetic light curve is short, ~ 13.5 
days. We compare this syn thetic light curve to the UVOIR 
light curve of SN2003hv (|Leloudas et al.1 120091 ) in Fig. 11, 
and find good agreement in shape from the earliest obser- 
vations, which are near maximum, to the latest ones, ~ 800 



days later. The rising part of the light curve is not con- 
strained by the data. Since we could derive the distribution 
of 56 Ni from the nebular spectra only, we have not included 
any 56 Ni at v > 10000 km s -1 . If some 56 Ni was present 
at higher velocities, the rising part of the light curve would 
be smoother, and the peak may be reached a little earlier. 
Figure 11 also shows a model where both the 56 Ni mass 
and the total mass located at velocities v < 8000 km s _1 
are multiplied by a factor 1.25, reflecting inaccuracies of the 
nebular-epoch modelling (see Section 5.1). The light curve 
for this model reaches peak somewhat later (~ 14 days after 
explosion) and is overall ~ 0.2 mag brighter. In both cases 
the inconsistency between the peak and the tail of the light 
curve is eliminated. 

The reason for the rapid rise is that not only the inner 
ejecta have low density, and that the total mass of Fe-group 
species is quite small (O.24M ), which keeps the opacity 
low. This is compounded by a normal ejecta velocity (the 
maximum-light spectrum of SN 2003hv has a Si II line ve- 
locity of ~ 10500 km s _1 ). Although what we have used here 
is by no means a realistic explosion model, the result that 
we obtained is nevertheless suggestive. 



6 CONCLUSIONS 

The high ionization and low near-IR flux in the nebular 
phase of SN 2003hv can be best explained if this SN la 
ejected less mass at low velocities than in standard single- 
degenerate, Chandrasekhar-mass explosion models. In par- 
ticular, the mass of stable NSE elements is greatly reduced. 
The smaller mass may lead to a short rise time, which would 
also solve the apparent inconsistency between the 56 Ni mass 
derived from the peak of the light curve and the nebular 
epoch. Therefore, the 56 Ni mass obtained from the peak lu- 
minosity using a relation which has been calibrated for a rise 
time of 19 days may be incorrect. 

Two alternative scenarios exist to the single-degenerate, 
Chandrasekhar-mass one. One is the explosion of sub- 
Chandrasekhar mass white dwarfs. This involves a CO white 
dwar f accreting He from a companion and an edge-lit explo- 
sion ilLivndll99Q],|Wooslev fc Weaverlll994l ; ILTvne fc Arnettl 
Il995l ; iFink et allboiOT ). ISim et all (|2010h calculated light 
curves of sub-Chandrasekhar models which do not consider 
how the explosion is triggered. Their more massive models 
predict risetimes of ~ 18 days. In those models, however, Fe- 
group elements are much more abundant than in the model 
we have derived, which is likely to result in a larger optical 
opacity. Their least massive model (O.88M ) has a risetime 
of 14 days, which is comparable to our result. In this model 
Ni is more concentrated towards the centre th an the dis- 
tribution we have derived for SN 2003hv (Fig. 7; ISim et al.1 
120101 . Fig. 1), and the kinetic energy is smaller (8.6 ■ 10 50 erg 
compared to ~ 1.25 • 10 51 erg). The diffusion time of optical 
photons may therefore be comparable in SN 2003hv despite 
the slightly larger mass. 

The other scenario is the mergin g of two white dwarfs 
l|lben fc Tutukovll9"83 ; IWebbinkll984l). and in particular the 
case of dynamical mergers jPakmor et aLll2010l ). This results 
in an explosion ejecting at least a Chandrasekhar mass, but 
the two merging white dwarfs had mass smaller than the 
Chandrasekhar mass, and consequently lower central densi- 
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Figure 12. A comparison of various density profiles: the W7 
density distribution (grey); the modified model with a reduced 
density in the inner layers obtained from our nebular calculations 
(red); a model of an explosion following t he dynamical merging 
of two white dwarfs llPakmor et al. || 201fj) f black), two different 
sub-Chandrasekhar models fblue) jsirn et alj|20ich . 

ties. This is reflected in the density distribution of the SN 
ejecta. 

In Fig. 12 we compare the density stru ctures of W7, 
of a merger model with total mass 1.8M© (|Pakmor et al.l 
l20ld ). of two sub-Chandrasekhar models, with total mass 
0.88 and 1 .06M W and 56 Ni mass of 0.07 and 0.43M o , re- 
spectively l|Sim et al.l l20ld ) , and the density structure we 
derived from our nebular model. We stress that our model 
is reliable only inside of 10000 km s -1 . In the inner ejecta 
(v < 10000 km s _1 ) the two sub-Chandrasekhar models and 
the density distribution we derived actually look quite simi- 
lar. A model with 56 Ni mass similar to that of SN2003hv is 
not available. 

It will be interesting to submit all non-standard models 
to the test of nebular spectroscopy. Clearly, the availability 
of flux-calibrated near-IR spectra at late times is an impor- 
tant piece of information if we want to discriminate among 
different models. 

The exact incidence of SNe with a non-standard origin 
is not known. In most SNe la, optical nebular spectra do not 
suggest a high Fe m/Fe 11 ratio. A dominance of Fe in over 
Fe n is deduced also in the innermost core of SN 1991bg, 
for which an improved fit to a series of nebular spectra can 
also be obtained by slightly reducing the mass at low ve- 
locities (Mazzali et al., in preparation). Both SNe 2003hv 
and 1991bg are characterised by a low luminosity, although 
the former is still spectroscopically normal at maximum. A 
possibility is that all SNe la at the dim end of the distribu- 
tion are the result of non-standard explosions, possibly of 
sub-Chandrasekhar mass white dwarfs. 

The final question is whether the unexpectedly low 
near-IR flux, which is one of the characteristics that leads to 
a low mass estimate for SN 2003hv, is typical of other SNe la. 
Late-time near-IR spectroscopy is available only for tw o 
other SNe la, 2003du and 2005W (|Motohara et al.l l2006h . 
There are also a couple of photomet ric observations in the 
near-IR, for SNe2000cx and 2001el jSollerman et al.ll2004 
IStritzinger fe Sollermanl l2007i ). In both of these cases, the 
ratio of near-IR to optical flux indicated a lower near-IR 
flux than expected from standard models like W7. While it 



is possible that other SNe may be accommodated in scenar- 
ios where the inner mass is reduced with respect to classi- 
cal models based on this evidence, we refrain from making 
that suggestion here. Near-IR and mid-IR spectrophotome- 
try would be required to be more confident about any such 
statement. In particular, SN 2000cx was very bright, indicat- 
ing that most of the inner part of the white dwarf progenitor 
had been burned to NSE. Also, spectroscopic modelling did 
not in dicate the need for a low density in the innermost 
layers jMazzali et alj|2007af ). 

In conclusion, near-IR spectrophotometry of SNe la 
in the nebular phase has proved to be a very important 
source of information in the case of SN 2003hv. It is to be 
hoped that such observations will be performed routinely for 
nearby SNe la (as well as other Type I SNe) in the future. 
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